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Abstract 

Light propagation in all-dielectric rod-type metamaterials is studied theoretically. The electric 
and magnetic dipole moments of the rods are derived analytically in the long-wavelength limit. The 
effective permittivity and permeability of a square lattice of rods are calculated by homogenizing 
the corresponding array of dipoles. The role of dipole resonances in the optical properties of the 
rod array is interpreted. This structure is found to exhibit a true left-handed behavior, confirming 
previous experiments [L. Peng et al, Phys. Rev. Lett. 98, 157403 (2007)]. A scaling analysis 
shows that this effect holds at optical frequencies and can be obtained by using rods made, for 
example, of silicon. 
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Metamaterials (MMs) are artificial structures made of microscopic elements whose col- 
lective behavior results in unusual macroscopic optical properties j^JLJsj]. Great efforts 
are currently being made to scale MMs down to optical frequencies Such an 

achievement would indeed make possible the development of novel, sophisticated, optical 
technologies in many areas, including telecommunications, life sciences and solar power. 

I 111 . [I2j have made a step forward by proposing to use high- 
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permittivity dielectric objects instead of metallic ones to avoid the losses and saturation 



effects inherent to the metal in the optical range [13|. This approach relies on the resonant 
modes that dielectric objects support [jjl]. Collections of resonators are expected to strongly 
modify the propagation of light at frequencies close to the resonances. Previous studies 



have for example noticed a correlation between the resonances of sing 
and the opening of photonic band gaps (PBGs) in arrays of them 15 



e dielectric objects 



161 ] . In the context 



of MMs, arrays of dielectric rods in p-polarized light (magnetic field parallel to the axis 
of the rods) have been shown to possess an effective, dispersive, magnetic permeability 
[17I . lEjj ]. Interestingly, it has further been suggested that dielectric rods in s-polarized light 
(electric field parallel to the axis of the rods) could exhibit both electric and magnetic 
dipole resonances, and thus, possibly constitute a so-called left-handed (LH) medium with 



simultaneously negative permittivity and permeability [111 . |l2j. These dipole resonances 
have been explained in terms of strong charge displacements and displacement currents, 
but to our knowledge, no theory on the electric and magnetic dipole activities of dielectric 
rods and their role in the optical properties of rod arrays has been given up to now. This, 
however, is an important matter if such structures are to be used in practical MM-based 
applications operating at optical frequencies. 

In this Letter, we present a theoretical study on the optical properties of periodic arrays 
of dielectric rods from the point of view of MMs. Our objective is to prove that these 
structures offer a similar control over light as conventional metallic MMs and that typical 
MM properties (e.g. left-handedness) can be observed at optical frequencies in very realistic 
and simple designs. We proceed as follows. First, we describe the electric and magnetic 
dipole activities of isolated rods by deriving explicit expressions of the corresponding dipole 
moments in the long wavelength limit. Second, we compute the effective material parameters 
of a square lattice of rods from these expressions and explain how the electric and magnetic 
dipole resonances of the rods cause the appearance of PBGs and LH dispersion curves. 
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Third, we show that these effects can be reproduced at different frequencies by tuning the 
rod resonances and illustrate our claim by the numerical demonstration of a true LH behavior 
at optical frequencies. 

We start by considering an isolated, infinitely long dielectric rod of circular cross-section 
C, radius R and relative permittivity e, surrounded by air, in a cartesian coordinate system 
of unit vectors u x , u y and u z , where the rod axis is along the latter. A planewave of 
wavevector k (|k| = k = 2ir/\) propagating along the x-direction illuminates the rod. The 
following theory is given for s-polarized light but similar steps could be carried out in the 
p-polarization. 

The scattering of light by circular cylinders is described by Mie theory, which provides 
exact analytical solutions of Maxwell's equations [l4]. In particular, the scattered electric 
field E s in the far zone (kr 3> 1) is given by: 

/FT ikr ( +°° \ 

-4j= e_<5 6o + 2^6 n cos(7i0) u, (1) 



7T v kr . 

\ n=l 

where 9 is the angle with respect to the x-direction and b n the nth-order Mie scattering 
coefficient of the rod. In fact, the polarization per unit volume P = Eq(e — 1) E induced 
in the dielectric rod acts as a source for the scattered field. E s can then be written in an 
integral form using Green's theorem as jl9]: 

„2 
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^ f c H^\k\v-r'\){e~l)W)d 2 r' (2) 



where Hq is the zeroth order of the Hankel function of the first kind. The integral has been 
restricted to C because the polarization vanishes outside the rod. 

A second expression of the far-field expansion in Eq. (CQ) can now be derived by expand- 
ing Eq. into a series of multipoles. Note that this technique differs from the familiar 
three-dimensional multipole expansion of the magnetic vector-potential [^(J because the 
bidimensionality of our problem necessarily implies a strong effect of the light polarization 
on the scattered field. In the far zone, Hq 1 ^ can be expressed by its asymptotic form 21]. 
The multipole expansion is introduced by writing |r — r'| ~ r — u r ■ r', where r = ru r , 
yielding the approximations \fk\r — r'| ~ \fhr and e* fc ' r_r '' ~ e lhr ■ e~ lhUr ' r ' . The expo- 
nential e _ifcUr ' r ' is then expanded in series of the source extension versus the wavelength as 
g-ifeur-r' _ ^~ tfc "^' r - ■ By inserting these expressions in Eq. (j2j), we obtain the polarized 
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multipole expansion of the scattered electric field in the far zone: 

[2 e ikr 

7T \fhr 



n=0 

with f n (r) = f t^L J c ( Ur ■ r')" (5 - l)E(r')dV. 

The successive terms of this expression describe 2D multipole radiation fields at large 
distances, the zeroth (n = 0) and first (n = 1) orders being assimilated to electric and 
magnetic dipoles, respectively. The corresponding scattering orders can be written as a 
function of the electric and magnetic dipole moments per unit length, defined respectively 
by p = j c P(r')d 2 r' and m = |/ c r' x J(r')d 2 r', with J = dP/dt the polarization current 
density. By equating them with their counterpart in Eq. (PQ), we can write p and m as a 
function of the Mie scattering coefficients b and & a , respectively, as: 

p/e = (4b /ik 2 )u z (4a) 

mZ = (-AbJik^Uy (4b) 



with Z = y/io/^o the free space impedance. These two expressions describe analytically 
the dipole activities of isolated rods in the long- wavelength limit in terms of their scattering 
matrix. 

The Mie scattering coefficients bo, b\ and 62 are plotted in Fig.[T]for rods with permittivity 
e = 600 llj. On the whole, |6q| and |&i| are found to remain larger by a few orders 
of magnitude than the higher-order coefficients, suggesting that the electric and magnetic 
dipole activities of the rods contribute to the most part of the optical properties of rod 
arrays. Higher-order multipoles, which are expected to induce spectrally narrow optical 
features near their resonance frequencies, can safely be ignored. 

We now consider an array of dielectric rods. In the long- wavelength limit, the corre- 
sponding array of dipoles can be described as an effective medium with permittivity and 
permeability dyadics E and p,, respectively. These macroscopic parameters are calculated 
from the microscopic polarizabilities of the rods through a process of homogenization, which 
takes the density of resonators and their mutual interaction into account. In the case of 
s-polarized light propagating along the x-direction, only the e zz and fi yy components are 
required to define the refractive index of the effective medium n e R = ^/E zz {i yy . Consider- 
ing that the incident electric field amplitude has been normalized to unity and using the 



4 




FIG. 1: (Color online) Complex modulus of bo (thick blue solid line), 61 (thin gray solid line) and 
62 (thin red dashed line) for rods with permittivity e = 600, in logarithmic scale. The insets show 
the amplitude of the scattered electric field at the first maxima of | bo | and | b\ \ . 

relations |H*| = |E*|/Zo and Eqs. (j4j), the electric and magnetic polarizabilities per unit 
length of the rods can be written as a e zz = p z /e El = Abo/ik 2 and a™ = m y /H l y = ibi/ik 2 , 
respectively. 

The effective material parameters of a square lattice of rods with permittivity e = 600 
and radius R = 0.68a/3, where a is the lattice per iodicity, are calculated using the nonlocal 
homogenization model proposed by Silveirinha [22J in the approximation of wavevectors close 
to the T-point. For the sake of comparison, this structure is similar to the one studied in 
Ref. As shown in Fig. [2] (left), the first electric dipole resonance of the rods [first inset of 
Fig. [1] induces a strong resonance of the permittivity, which results in the opening of a wide 
frequency range of negative permittivity. At a/A ~ 0.07, the first magnetic dipole resonance 
of the rods [second inset of Fig. [1] induces a sharp resonance of the permeability, which 
eventually takes negative values within the neg ative permittivity range, thereby inferring a 



LH behavior in accordance with Peng et al. 



In Fig. [2] (right), we compare the dispersion 



curves of the effective medium calculated using the relation q x = n c ^uo/c and those of the 
rod array, calculated by the planewave expansion (PWE) method using a freely available 



software package 23J. Apart from symmetry degeneracies and higher-order resonances, 
which have not been considered in our theory, the main optical features of the rod array are 
very well reproduced. We can then confidently say that this structure truly behaves as a 
medium with negative permittivity in the two lower-frequency PBGs, and as a LH medium 
in the narrow frequency range at a/A ~ 0.07. 

At this point, we have shown that the optical properties of periodic arrays of rods result 
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FIG. 2: (Color online) (Left) Real parts of the effective permittivity e zz (thick blue line) and 
permeability \i yy (thin gray line) of the rod array (R = 0.68a/3, e = 600). (Right) Dispersion 
curves of the structure calculated from the effective material parameters (black solid lines) and by 
the PWE method (gray dashed lines). 

from the collective response of the resonant rods. These structures rely on similar principles 
as conventional metallic MMs and thus, may be used as such to control the propagation of 
light. 

Due to the growing interest in developing MMs for the optical range, it is now important to 
investigate the scaling properties of all-dielectric rod-type structures. Previous studies have 
limited their work to high-permittivity rods to place their resonances in the homogeneous 



3 nomogen 
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regime (A 3> a, R) and prevent them from exhibiting a strong spatial dispersion [8|, |ll|, [12]. 
Our theory provides additional information on this matter. As shown above, the electric and 
magnetic dipole activities of dielectric rods are intrinsically related to their Mie scattering 
coefficients. The scaling properties of rod-type structures may therefore be understood by 
studying the variation of these coefficients with the rod permittivity e. Decreasing e reduces 
the optical size of the rods and thus, the wavelength A at which they resonate. This is 
evidenced in Fig. [3] (left), where the wavelength positions of the electric and magnetic dipole 
resonances (indicated by the maxima of |&o| and are shown to depend almost linearly on 
the rod refractive index n = ^fe. In particular, the magnetic dipole resonance observed at 
X/R ~ 63 (a/A ~ 0.07) in rods with permittivity e = 600 is shifted to X/R ~ 8.8 (a/A ~ 0.5) 
when e = 12. As shown in Fig. [3] (right), this permittivity is sufficiently high for 6 and b\ to 
remain the most significant coefficients up to the resonance frequency of b 2 at Rj A ~ 0.17. 
The LH behavior is therefore expected to hold. 

This is verified by comparing the photonic band structures and second-band iso-frequency 
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FIG. 3: (Color online) (Left) Wavelength positions (in units of X/R) of the two lower-frequency 
maxima of |&o| (thick blue lines) and \bi\ (thin gray lines) as a function of the rod refractive index 
n. (Right) Complex modulus of bo (thick blue solid line), b\ (thin gray solid line) and 62 (thin red 
dashed line) for rods with permittivity e = 12. 

curves (IFCs) of the higher and lower-e structures calculated by the PWE method. As 
observed in Figs. H^a)-(d), their optical features are readily similar. First, the LH dispersion 
curves lying at a/A ~ 0.07 are pushed up to a/A ~ 0.5, as expected. The splitting between 
the two higher-frequency bands is larger in the latter structure because the decrease of the 
rod refractive index broadens the rod resonances [Fig. [3] (right)] and increases the inter- 
rod interaction [15] . Second, the IFCs both exhibit a strong spatial dispersion even in the 
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251 ]. arguing that 



vicinity of the T-point. This observation supports previous studies 
large wavelength-to-period ratios do not necessarily result in isotropic responses. The spatial 
dispersion here is naturally inferred by the magnetic dipole activity of the structure H] ■ 
The lattice symmetry therefore plays an important role in the spatial response of the rod 
array in the LH frequency range. 

To evidence the existence of a LH behavior at optical frequencies, we finally simulate 
the propagation of light in the lower-e structure with the 2D finite-difference time-domain 
method using a freely available software package 



Light is incident at an angle of 20 c 



and the structure is tuned to near- infrared wavelengths (A = 1.55 /im) by using rods of 
radius R = 158 nm and a lattice of periodicity a = 698 nm. As shown in Fig. 11(e), the 
phase of the propagating field in the rod array is opposite to that of the field in free space, 
which indicates, as expected, a LH behavior. 

Taking all the above into consideration, we come to the conclusion that the LH behavior 
initially observed by Peng et al. [11[ is not specific to large wavelength-to-period ratios, for 
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FIG. 4: (Color online) Photonic band structures (a, c) and second-band IFCs (b, d) of the higher 
and lower-e rod arrays (top: e = 600, down: e = 12). The IFCs are labeled by their reduced 
frequency a/A. (e) Steady-state amplitude of the electric field of light at A = 1.55 /jm incident at 
20° on the lower-e structure (i? = 158 nm, a = 698 nm). 

this effect is primarily a matter of coupled resonances. 

It is particularly important to note that the physics involved here remains the same 
even though, in lower-e structures, A becomes comparable to a. In this regard, it would be 
interesting in future studies to investigate to what extent these structures can be considered 
as homogeneous. Concurrently, our theoretical results could be used to clarify the origin of 
the many LH behaviors that have been reported in rod-type photonic crystals in the past 



few years (see e.g. Ref. [291] ). Experiments could also be initiated at once, as techniques 
to fabricate silicon-based (e ~ 12) rod-type structures already exist |3(|. Silicon holds a 
pre-eminent position in Photonics and is therefore a very interesting candidate material for 
all-dielectric MMs operating at optical frequencies. Rod-type structures are finally likely 
to integrate well on photonic platforms, considering that high couplin g ef ficiencies can be 



obtained even at large angles of incidence owing to the rod resonances 



3l|. 



To conclude, we have presented a theoretical study on the optical properties of all- 



S 



dielectric rod-type structures, using an approach based on the electric and magnetic dipole 
activities of dielectric rods. This work constitutes a first proof that dielectric rods can be 
used to design true MMs at optical frequencies. Due to the large experimental knowledge 
in fabricating nanoscale dielectric structures, exciting applications, such as all-dielectric in- 
visibility cloaks, may be realized in the optical domain in the near future. 
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